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ABSTRACT 

The microlensing probability (optical depth r) toward the Galactic center 
carries information about the mass distribution of the Galactic bulge/bar, so 
can be used to constrain the very uncertain shape parameters of the Galactic 
bar. In this paper we study a family of bar models with tunable triaxiality, 
boxyness and radial profile to span the range of plausible density models of 
the Galactic bar. It also includes four models that were used to fit to the 
COBE maps by Dwek et al (1995). We find the optical depth on the minor 
axis of the whole family of models and in special cases the whole microlensing 
map of the bulge/bar are given by simple analytical formulae. The formulae 
show the following dependence of the optical depth on the bar mass, radial 
profile, angle, axis scale lengths and boxyness. 

(1) t is proportional to the mass of the bar, M. 

(2) r falls along the minor axis with a gradient determined by the density 
profile and boxyness of the bar. 

(3) Going from an oblate bulge to a triaxial bar, the optical depth in- 
creases only if a < 45°, where a is the angle between the bar's major axis 
and our line of sight to the center. 

(4) Among bars with different degree of triaxiality, r is the largest when 
the bar axis ratio |p = tana, where yo,xo are the bar scale lengthes on the 
short and long axes in the Galactic plane respectively. 

(5) At a fixed field on the minor axis but away from the center, boxy bars 
with a flat density profile tend to give a larger optical depth than ellipsoidal 
bars with a steep profile. 

(6) Main sequence sources should have a significantly lower (20-50% lower) 
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optical depth than red clump giants if main sequence stars are not observed 
as deep as the bright clump giants. 

For the four COBE-constrained models, the optical depths contributed 
by the bar is in the range of w (0.3 - 2.2) x lCT 6 M/(2 x 10 10 M Q ) at Baade 
window. Out of the four models, only the model with a Gaussian profile and 
ellipsoidal shape satisfies the observational constraint, the other three models 
produce optical depths which are inconsistent or marginally consistent with 
the 2(T lower limit of the observed optical depths, even if we adopt both a 
massive bar 2.8 x 1O 1O M0 and a full disk. Independent of the COBE map, 
we find that the microlensing models can potentially provide constraints on 
the mass distribution. If increasing microlensing statistics confirm the high 
optical depth, ps 3 x 10~ 6 , presently observed by the MACHO and OGLE 
collaborations, then the observation argues for a massive (> 2 x 10 10 M Q ) 
boxy bar with axis ratio ^ ps tan(a) and a < 20° and with a flat radial 
profile up to corotation. Stronger limits on the bar parameters can be derived 
within two years when the sample is increased from the current fifty events 
to about two hundred events. 

Key words: dark matter - gravitational lensing - galactic centre 

1 INTRODUCTION 

Microlensing surveys such as the MACHO (Alcock et al. 1993, 1995a, b), OGLE (Udalski 
et al. 1993, 1994), EROS (Aubourg et al. 1993) and DUO (Alard et al. 1995) collaborations 
have discovered more than 100 microlensing events toward the Galactic bulge and about 
10 events toward the LMC. One of the many exciting discoveries (see Paczyhski 1996 for a 
review) is the surprisingly high optical depth, r, i.e., microlensing probability, toward the 
Galactic centre (Udalski et al. 1994; Alcock et al. 1995b). The OGLE collaboration found 
r_ 6 = t/10~ 6 = (3.3 ± 2.4) (2a error bar) based on 9 events and the MACHO collaboration 
found t_ 6 = 2.43i ; 74 (2a) based on the full sample of 41 events and r_ 6 = 6.32l 3 57 (2a) 
for 10 low latitude clump giants. This is clearly in excess of all the predictions prior to the 
discovery of microlensing events (Kiraga & Paczyhski 1994a). Both a maximum disk (Gould 
1994a, Alcock et al. 1995b) and a Galactic bar oriented roughly toward us (Paczyhski et al. 
1994a; Zhao, Spergel & Rich 1995, Zhao, Rich & Spergel 1996, hereafter ZSR95, ZRS96) 
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have been proposed as an explanation. The latter interpretation is favored by other lines of 
evidence, such as the non-circular motion of gas in the inner Galaxy (Binney et al. 1991; 
Blitz & Spergel 1991), the asymmetric and boxy COBE map (Weiland et al. 1995) and that 
the clump giants appear brighter at one side (/ > 0°) than the other side (7 < 0°) of the bulge 
(Stanek et al. 1994). A Galactic bar, together with a normal disk, can produce an optical 
depth r_6 ~ 2.5 (ZSR95, ZRS96). This value is in good agreement with the OGLE and the 
MACHO full sample values, but seems to be at the la lower limit of that for the 10 low 
latitude clump giants in the MACHO sample. The theoretical prediction depends on the still 
uncertain bar axial ratio, orientation and mass profile; Dwek et al. (1995) found that several 
bar volume density models are consistent with the COBE map. Given these uncertainties, 
what is the range of the predicted optical depth? What are the optimal bar configurations 
that can produce an optical depth as high as r_ 6 = 6? Can the Galactic microlensing 
observations put constraints on the bar parameters? These are the main questions we want 
to address in the paper. We demonstrate that the published optical depth based on about 
10 microlensing events with red clump giants as sources already put interesting constraints 
on the orientation, axial ratio and boxyness of the bar. With steadily increasing number 
of events and decreasing error bar for the optical depth, microlensing will become a very 
powerful tool to probe the Galactic structure. 

The outline of the paper is as follows. In §2, we study a family of models for the bar 
density distribution, which includes many models used in the literature. We derive an ana- 
lytical expression of their optical depths on the minor axis and for ellipsoidal bars the whole 
microlensing map on the sky as well. In §3, we apply our formalism to the Galactic bar, 
and compare the results with the observations. We summarize our results and discuss the 
implications in §4. More mathematical details are given in the appendices and at the end of 
Appendix C some instructions of using our formulae. 

2 MODEL 

In this section we first give the general expressions for the microlensing probability (optical 
depth) of a bar, then derive the optical depth map for a simple ellipsoidal Gaussian bar 
model. The results are then extended to a family of bar models with more general shapes 
and radial profiles. Using these models we study the dependence of the optical depth on 
parameters of the bar. 
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2.1 Microlensing Optical Depth of the Bar 

The microlensing optical depth for a source at distance D s is the probability for it to come 
within one Einstein radius of any lens (deflector) placed at < Dd < D s : 

r ( D.) -J* (^>) (,#,) dD d , RE .(^Ef, fls(fl ,- fl ^, (1) 

where m is the mass of the lens, p(Dd) is the mass density of the lens at distance Dd, Re 
is the Einstein radius, D is the effective distance between the lens and the source, c is the 
speed of light and G is the gravitational constant. The first bracketed term is the number 
density of the source at distance Dd, while the second term gives the lensing cross section, 
namely, the area enclosed by one Einstein radius. Note that the lens mass m cancels out, 
therefore the optical depth depends only on the total mass density of the lenses, but not on 
the mass spectrum. 

The microlensing optical depth is usually measured for all sources in a narrow cone of 
the sky (e.g., Baade window of the bulge, / = 1°, b = —3.9°), therefore it is necessary to take 
into account the distribution of source distances by averaging over the line of sight: 

where D s max is the maximum distance of an observable source, and n(D s )dD s is the number 
of observable sources at distance D s to D s + dD s in the field of view. The factor of D 2 S in 
n(D s ) takes into account the volume increase with distance for a fixed field of view, and 
the factor D~ 2f3s parametrizes the decrease of detectable sources with distance due to the 
detection threshold (Kiraga & Paczyhski 1994). For the main sequence stars, Kiraga & 
Paczyhski (1994) and ZSR95 suggested 1 < /3 S < 2. For the clump giants, (3 S = is more 
approriate as they are probably seen throughout the Galaxy. Note that (r) is generally a 
function of Galactic coordinates, but we ignore the angular variation inside a small field of 
view (less than one square degrees). 

In all the following sections, we will use the Galacto-centric coordinate system (x,y,z), 
with the Galactic centre at the origin, and the solar system at (—R ,0,0), where Rq = 8.5 
kpc is the distance to the Galactic centre. (X, Y, Z) are used to denote the coordinates along 
the three principal axes of the bar; X being the major axis of the bar. Since the tilt out of 
the Galactic plane is negligible (Dwek et al. 1995), the (X, Y, Z) coordinate system is simply 
related to (x, y, z) by a rotation of angle a in the x — y plane. It is understood that only the 
absolute value of the angle, \a\, is used. 
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Since the longest dimension of the bar, a, is still much smaller than the bar's distance 
from us, Ro, we can make two simplifications to eq. ([|) so that the calculations are tractable 
analytically: 

(1) We set n(D s ) = p(D s ), i.e., j3 s = 1, so that the volume increase with distance just 
balances the decrease of detectable sources. This is a plausible assumption for the bulge 
main sequence stars, but probably underestimates the optical depth for the clump stars 
with f3 s = by an order of (■§-)• 

(2) We adopt a plane parallel approximation: 

x s & D s - R , xi^D d -R , ^ « 1, D = (D s - D d )^- « x s - x h 
and the Galactic longitude and latitude of any bulge star 

M)-(f,f). 

Here (y, z) are the coordinates where the line of sight intersects with y — z plane at the 
center. x s and xi denote the ^-coordinates of the source and the lens, which are assumed to 
be in the range of 

— OO < X[ < x s < +oo. 

In making this simplification, we effectively set — Rq — > — oo, and s max — > oo. This is 
valid for lenses in the bar as long as the density of the bar falls off sufficiently fast or is 
truncated within a few kpc, which is the case for all our models. Note by letting = 1 
we overestimate D, and hence the optical depth, by an order of (^)- The assumption also 
neglects the difference in distance between the positive and negative longitude sides. 

With these simplifications, the optical depth for a source located at (x s , y, z) is given by 
(cf. eq. 1) 

AtiG f Xs 

t(x s , y, z) = — r- / p{x h y, z) (x s - x t ) dx h (3) 

C J — oo 

while the average optical depth toward a field (y, z) is (cf. eq. 2) 

/ / z ^ = 1-^ dx s p{x s , y, z) t{x s , y, z) 
I-^dx s p(x s ,y,z) 

For the clump giants, the effects of our assumptions (1) and (2) on r are in the opposite 
direction. Remarkably, numerical calculations find that these two effects cancel out, there- 
fore, our analytical results based on eqs. ([|) and (Q) are accurate for the clump giants (cf. 
Fig. 3.) For the fainter main sequence stars eq. ([|) overestimates the optical depth. But this 
can be easily corrected analytically to the first order of using eqs. ( |C21|) and (|C20| ) in 
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Appendix C. The corrected values are again in good agreement with numerical results. The 
analytical formulae, of course, have the significant advantage that they allow the dependence 
on the parameters to be studied explicitly (see §2.3.) 

2.2 Analytical Results For a Gaussian Ellipsoidal Bar 

Consider the optical depth of the Galactic bar with a Gaussian radial profile and an ellip- 
soidal shape, which is one of the models that Dwek et al. (1995) used to fit the COBE map. 
The density of this model is given by 

pO, V, z) = po exp l-^r) , Po = - .3 M , (5) 

V (27c)2x y Zo 

where 



A \x Q ) + \y ) + \Zq) ' ^ 
M is the total mass of the bar, xo,yo, zq are the scale lengths along the three principal 
(X, Y, Z) axes of the bar and po is the central density. X, Y, Z are related to x, y, z by a 
rotation of the bar angle a, 

X = x cos a — ysina, Y = xsina + ycosa, Z = z. (7) 
Substituting eq. (|7|) into eq. (|6|), one finds that 



A 



2 



X 



where Xi, y\ and Z\ are effectively the scale lengthes along the x, y and z axes. Together 
with k they are given by 

1 cos 2 a sin 2 a 1 sin 2 a cos 2 all 11 

O o H - o i o 9 ~~ \~ o 5 9" 9~ 5 ^ \ 9 9 J Sill Ot COS Ol . ( C/ ) 

x i x o yo yi x o yo z i z o yo x o 

So for the Gaussian ellipsoidal bar, the optical depth for a source at the Galactic centre 

is (cf. eq. ^) 

AnG ro . x? . . , , AttG 

If y_ oo Poexp(-^-2)(-xO^ = — 

The optical depth for a source on the minor axis at (0, 0, z) is given by 

r(0,0,^)=r(0,0,0)exp(-^). (11) 

iz 

The optical depth for all sources along a line of sight passing the minor axis can be obtained 
from eq. (|J) 

i- , x 4:irG M 1 , yn s / z 2 . , . 

(r(0, z)) = V2t(0, 0, z) = — --^lA a) exp(-- ? ), (12) 

© 1996 RAS, MNRAS 000, 



r(0, 0, 0) = —j- / Poexp(-— = — ^-poz?- (10) 



where 

T ,yo x 

J(— ,a) 
x 



1 

2^ 



yo 2 , 

— cos a ) + 

x 



Xo . 2 

— sin a 



-i 



1 (u + u- 1 )' 

7T 



u 



yo/x 



(13) 



sin 2a tan a 

Similarly we can obtain a microlensing optical depth map on the whole sky, i.e., on the 

whole y — z plane: 

( r (y, z )) 



A-rGM 1 y 

-j^I{—,a)exp 



C 2 Zq-K^^X^'^ 2yl 2 



2 Z r 



x\ sin 2 a + yl cos 2 a ( 14) 



Note that for any field with the same z, the optical depth is the largest on the minor axis. 
2.3 Analytical Results For a Family of Bar Models 

The bar density is not uniquely constrained by the COBE map (Dwek et al. 1995). Models 
with various radial profiles and shapes can fit the data about equally well. In this section, 
we derive the optical depth for a general set of bar density models and show its dependence 
on the profile and shape of the bar. 

Consider the following class of bar density models: 

p(x, y, z) = p f(X), /(A) = fl + —j , (15) 

where the central density po can be expressed in terms of the total mass of the bar, M, and 
is given in eq. (Al) in the Appendix A. The effective dimensionless radius parameter A is 
given by 



A' 
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k 

+ 


Y 


k 


k 

+ 


z 




X 




yo 







(16) 



The parameters n and k determine the shape of the bar: k specifies the shape of the bar in 
the X — Y plane while n specifies the shape in the Y — Z and X — Z planes. The parameter 
13 describes the radial profile of the bar. When f3 > 0, the density falls off like a power-law 
with slope (— (3) at large radius. When (3 < 0, the model describes a truncated finite bar 
confined within A < \/3 — 3|«. When ft — > oo, one obtains an exponential or Gaussian density 
profile for n = 1 or n = 2. The Gaussian ellipsoidal model studied in the previous section 
corresponds to a special case with (/3,n,k) = (oo,2,2). A homogeneous rectangular bar 
corresponds to a model with (/?, n, k) = (0, oo, oo); analytical results for this special case are 
given in the Appendix B. 

This family embodies the more interesting cases of bar models[j] and can simulate a great 

t models with < /3 < 3 correspond to flaring bars with their densities increase with radius, and will not be considered 
further. 
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diversity of shapes and radial profiles. The range of shapes and profiles are demonstrated 
by a few examples in Fig. la and lb. 

Similar to the Gaussian ellipsoidal bar model, one can show that for this family the 
optical depth at any minor axis field depends on z as: 

1 (\z 



<t(0,s)) 
where 



H(k, — , a) 
x 



AuGM y 
c 2 z x 



1 + 



(3-3 \zq 



K(k) 



Vo 2 \ = . (xq . 2 

— cos a) + I — sm a 

x ) \y 



fc-, -As/2 

2 



2 ( 



k 

U a 



li a 



-2/fc 



K(k) 



sin 2a 



(17) 



(18) 



£o((3,n) « 0.6 ± 0.15 is a slow varying function of /3 and n, u is as defined in eq. (|T3|). 
K(k) = 4k~ 1 B(k~ 1 , /c _1 ), and B is the complete Beta-function. For clarity, we have deferred 
the somewhat lengthy definition of £o (together with several useful approximations) until 
Appendix A. Eq. ( |T7D reduces to eq. (|T2D for a Gaussian ellipsoidal bar when (/3, n, fc) = 
(oo,2,2), with = ■K~ 1/2 ,H(2,y G /x ,a) = I(y /x ,a). 

For ellipsoidal models with n = k = 2, an analogous microlensing map as eq. fll4j) can 
be derived. We find 

1 

where is given by eqs. (|I3]) and 
is the largest on the minor axis y = 0. 

Eq. flTTD is our main analytical result. In the rest of this section, we will study the 



<T(y,z)) = <T(0,0)> 



1 + 



' I/ 2 ^ 2 



, ifn = Jfe = 2, (19) 
Along the line with a constant z, the optical depth 



dependence of the optical depth on the parameters. 



2. 3. 1 Dependence on Radial Density Profile and Latitude 

(1) The optical depth due to the bar lenses, r, is proportional to the total bar mass M, as 
expected. 

(2) t falls off on the minor axis as 



P(£) 
P(0) 



1-4 



9 (cf. eqs. [16| and |T7|). When j3 > 3, the model 
is a power-law with a finite mass, r falls off slower than the density. When (3 < 0, the model 
is truncated at a finite size, r falls off faster than the density. When (3 — >• oo, r falls off as 
the density. 

(3) For a fixed field z kpc above the plane, r first increases and then decreases with the 
vertical scale length zo, and r is maximized when the scale height zq = z. This is because 
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Figure 1. (a) Some examples of geometrical shapes are shown for the family of bar models as described by eqs. jl5[ ) and jl^). 
For illustration we used a bar with axis scale lengths xq : yo : zq = 2 : 0.8 : 0.5 (kpc), similar to the scale lengths of the COBE 
constrained GE model (cf. Table 1.) The shapes in the X-Y plane are determined by k while the shapes in the X-Z (also Y— Z 
but not shown) plane are determined by n. The shape of the bar evolves from being diamond, ellipsoidal, boxy to rectangular 
when k or n are increased in the sequence of 1, 2, 4, and oo. The dashed line indicates our line of sight to the center, which 
makes an angle a with the major (X) axis of the bar. Observationally, a ~ 10° — 45° with the near end of the bar in the first 
Galactic quadrant. 
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(b) Some examples of radial profiles are shown for the same family of bar models. The density profiles are described by two 
parameters n and /3 in the radial coordinates expressed by A (cf. eqs. ^ and All the density profiles are normalized to the 
central density. 
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increasing the scale height leads to both a shallower density gradient and a smaller central 
density. To obtain the maximum optical depth, we note that in eq. flT7|) the factor 



i+ 1 



(3-3\z 



-2 g-2 



M 

The equality sign indicates the maximum and is reached when z = z irrespective of j3. 



2.3.2 Dependence on Bar Angle and Shape 

Not surprisingly, r depends on some combinations of the angle, axis ratio and the shape of 
the bar, more specifically: 

(4) For fixed axis ratio and shape, r is largest when the bar points toward us, i.e., when 
a = 0. In this case, r oc H(k, 0) oc ^ (cf. eqs. and [1SD . 



(5) For a fixed angle a < 45°, gradually squashing an oblate bulge to make it increasingly 
elliptical in the x — y plane leads to a monotonic increase in the optical depth until ^ < 
tan(a). For a > 45°, the optical depth always decreases when Uo/xq changes from one to 
zero, r has a maximum at u — 1 (cf. eq. 18), i.e., at ^ = tana. This is graphically shown 



in the upper panel of Fig. 2 for k = 2 models. The turn-over for very small axis ratio can 
be intuitively understood as a needle-shaped bar pointing sideways has zero depth to place 
the lens and the source along the line of sight, hence zero optical depth. Mathematically, r 
is proportional to H (cf. eq. |l^), and 

„/, Vo v 2 2 2-i 

1 'x #(jfe) sin 2a — i^(/c) sin 2a' 1 ' 

So r and H are maximized at u = 1 irrespective of k. 

(6) The maximum gain in optical depth of an elliptical (k = 2) bar vs an oblate bulge is 
given by 

< -J— , (22) 
T bu i ge sin 2a 

where the maximum is at — = tana. For example a bar with a = arctan -1 — = 15° gains 
by a factor of 2. But for a > 45° the oblate model predicts higher optical depth than bar 
models. 

(7) Other things being equal, r is a slow varying function of n. The ratios for n — 1, 2, 4 
(diamond, ellipsoidal, boxy) for /3 — > +oo are about 1.3 : 1 : 0.95, while for j3 = 4 the ratios 
are 1.68 : 1 : 0.82. 

(8) Similarly, for fixed axis ratio and angle, varying k, r changes by less than a factor of 
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2 (see Fig. 2, lower panel). The ratios of the optical depth r for the diamond shape (k = 1) 
to that for an ellipsoidal shape (k = 2) are 1.57, 1, and 0.79 at u = ^(a) = 0, 0.31, and 1 
respectively. The variation is the same for u > 1 due to the symmetry between u and u _1 
(cf. eq. 0). 

(9) Of the whole class of models with different central densities and shapes, the optical 
depth is the largest for a homogeneous rectangular (HR) bar. The maximum optical depth is 
reached when a = tan" 1 (f^), i-e., when the line of sight coincides with the diagonal line of 
the homogeneous rectangular bar. The microlensing map of this model is given by equations 
in Appendix B. A comparison of the homogeneous rectangular (HR) model with models of 



other shapes is in § 3.5 



3 APPLICATION TO THE GALACTIC BAR 

The formalism we have developed can be applied to any edge-on bar. In this section, we 
will concentrate on the Galactic bar, where most microlensing events are observed (Alcock 
1995b; Udalski et al. 1994; Alard et al. 1995). As the optical depth is a function of the bar 
parameters such as the total mass, bar angle and axial scale lengths, we therefore first outline 
the uncertainties in these parameters before we proceed to estimate the optical depth. 

3.1 Uncertainties in the Bar Parameters 

Dwek et al. (1995) carried out a systematic study of the bar parameters by fitting the COBE 
map. A few models can fit the data well, i.e., the parameters are not well constrained. For 
non-axisymmetric models, the bar angle lies 10° ^ a ^ 45°, and the axis ratio Do/xq = 
0.15 — 0.5, and the vertical scale height z = 200 — 500 pc. In Table I, we have listed four 
models from Dwek et al. (1995) that can be described by our parametrization in eqs. ([15]) 
and flT6|) . Except for the power-law exponential model, the errorbars for the parameters are 
quite small (cf. Table 1), while the variations from model to model are quite large. For each 
model, we will take the best-fit (central) parameter values in the table and evaluate the 
optical depth. We do not attempt to do a full error analysis due to the correlation of errors 
between parameters. Out of the four models, the Gaussian ellipsoidal model provides the best 
fit to the COBE map, while the exponential diamond (ED), oblate bulge (OB) and power-law 
ellipsoidal (PE) profiles become progressively less satisfactory. Note that the oblate bulge 
model is an axisymmetric model. Two additional models are listed in Table 1, one being 
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Figure 2. shows the ratio of optical depths of a bar vs an oblate bulge as a function of the bar angle and shape. The upper 
panel shows that for an elliptical bar (k = 2), the optical depth peaks at an axis ratio = tan(o) if the bar angle a < 45°. 
For a > 45°, the optical depths increases monotonically with the axis ratio yo/xo, and is always smaller than that of an oblate 
bulge model. Note an elliptical bar becomes an oblate bulge at = 1. For a bar angle a = 15°, the lower panel shows that 
an elliptical bar gives a larger peak (maximum) optical depth than that for a diamond shaped bar (k = 1), but less than those 
for a boxy bar (k = 4), and a homogeneous bar. The positions of the maximum optical depths are indicated by triangles. They 
always occur at the axis ratio ^2. = tan(a) = 0.27 for a = 15°. 

a truncated boxy (TB) bar model, and the last one a homogeneous rectangular (HR) bar. 
These models are interesting since extragalactic bars often have a flat density profile and a 
sharp truncation near corotation (Sellwood & Wilkinson 1993). For these two models, the 
bar parameters are not available and we only consider their optimal lensing configurations. 
These six models should offer a broad survey of the uncertainty in the optical depth. 

The mass of the Galactic bar cannot be derived directly from the COBE map. However, 
it can be obtained by fitting either the gas rotation curve or the stellar velocity dispersions. 
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Table 1. Parameters for Six Galactic Bar Models 



Model 


OB (GO) 


GE (Gl) 


ED (El) 


PE (P3) 


TB 


HR 


Description 




Gaussian, 


Exponential, 


Power-Law, 


Truncated, 


Homogeneous 




Oblate 


Ellipsoidal 


Diamond 


Ellipsoidal 


Boxy 


Rectangular 


{/3, n, k) 


(oo,2,2) 


(oo,2,2) 


(oo,l,l) 


(4,2,2) 


(-4,4,2) 


(0, oo, oo) 


o 


— 


11 ± 2.4 


24.4 ±6.7 


45.4 ±33.3 


— 


— 


xq (kpc) 


0.91 ±0.01 


2.08 ±0.06 


1.64 ±0.06 


0.90 ±0.2 








Vo (kpc) 


0.91 ±0.01 


0.75 ±0.01 


0.61 ±0.11 


0.23 ±0.98 


— 


— 


zo (kpc) 


0.51 ±0.01 


0.45 ±0.02 


0.31 ±0.01 


0.28 ±0.04 






to{/3,n) 


0.564 


0.564 


0.75 


0.477 


0.554 


0.667 


U{P,n) 


1.772 


1.772 


2.55 


2.094 


1.376 


1 


H(k,^,a) 


0.159 


0.354 


0.165 


0.08 






(r-e(0,0)) 


2.9 


7.3 


6.5 


2.2 






(r_ 6 (0,-2°)> 


2.4 


5.8 


2.5 


1.1 






(r_ 6 (0,-4°)> 


1.5 


.3.1 


0.94 


0.41 






(t_ 6 (0,-6°)> 


0.64 


1.1 


0.38 


0.20 







Six models for a very massive (M = 2.8 X 1O 1O M0) Galactic bar. The radial profiles and geometrical shapes of the models are 
described in the second and third rows. The initials of these descriptions in the first row are used to label the curves in the 
figures. The first four models are the same as the GO, Gl, El, and P3 models in Dwek et al. (1995). (f3,n,k) are parameters 
that describe the shape of the bar (cf. eq. |l(j). a is the angle between the direction to the Galactic centre and the major (X) 
axis of the bar. xo, yo, zo are the scale lengths of the three principal axes of the bar. The parameters for the first four models 
are taken from Table 1 in Dwek et al. (1995) at 2.2 /im with R max = oo. §o(/3,n) , §4(/3,n) and H(k,yo/xo,a) are functions of 
f3, n, k, a and yo/xg. These are used in eqs. ( |l7j ) and (320) to evaluate the optical depth. The last four rows list the predicted 



optical depths due to lenses in the bar at four fields; the values would scale linearly for other choice of the bar mass M (eq. H 
gives the range of M). 

Generally axisymmetric models such as the Kent's (1992) oblate rotator model (cf. the oblate 
model in Table 1) fail to match the gas rotation curve either at R ~ 0.5 kpc or at R ~ 3 
kpc. Triaxial models with and without pattern rotation can be made consistent with the gas 
rotation curve. A rough limit on the mass of the bar can be set as follows: 

If the bar is contained within a corotation radius of R cor ~ 2.4 kpc (Binney et al. 1991), 
then the enclosed dynamical mass, Md + M, satisfies 

= V?(iW, (23, 

-ttcor 

where Md is the mass contributed by the disk within R COT , and V c (R COT ) ~ 195 km s _1 is 
the circular velocity at the corotation radius. The parameter r\ is a factor that describes the 
enhanced gravity of a flattened or barred system relative to a spherical system. Since t] > 1 
and Md ~ (10% — 20%)M > 0, we derive a loose upper bound for the bar mass, 

M ^ 2.2 x 1O 1O M (24) 

by setting r\ = 1 and Md = 0. 

More detailed gas modelling that matches the non-self-crossing closed orbits with the 
gas terminal velocity gives a very similar value. Assuming different pattern speeds, viewing 
angles, and bar density profiles, several authors find M = (2.2 ± 0.2) x 1O 1O M (Gerhard &; 
Vietri 1986; Binney et al. 1991; ZRS96) 

The mass obtained by different authors using the stellar velocity dispersions of bulge 
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stars has a larger scatter. This is mainly due to different levels of sophistication in the 
dynamical models. The simple oblate rotator model by Kent (1991) predicts a low value, 
M — 1.8 x 1O 1O M based on the pre-COBE data. The detailed steady state bar model of 
Zhao (1996), which fits the COBE map with a positive definite distribution function, gives 
M = (2.2 ±0.2) x 10 10 M Q . Straightforward but less certain application of the Virial theorem 
predicts between M = 1.8 x 1O 1O M with no pattern rotation of the bar (Han & Gould 1995), 
and M = 2.8 x 1O 1O M with a big pattern rotation (Blum 1995). 

To summarize, the predicted mass of the bar by different authors is in the range of 

1.8 x 1O 1O M < M < 2.8 x 1O 1O M . (25) 

3.2 Optical Depth of Disk Lenses 

To compare our bar models with observations, the contribution from the disk lenses must be 
added. It is well known that a standard disk model predicts Tdi S k,-6 ~ 0.66 at Baade window, 
too small to explain the observed high optical depth. However, our predictions for disk lenses 
are still uncertain due to lack of constraints of faint disk sources at several kpc from the 
Galactic center. ZSR95 find that the optical depth of the disk is in the range of 0.37 — 0.87 
for Baade window by extrapolating the local disk density with a double-exponential disk 
model with different disk scale length, scale height and local density normalization. The 
lower end is from a disk with a central hole as suggested by Paczyhski et al (1994b) based 
on color-magnitude diagrams of OGLE stars. A much higher disk optical depth (rdisk,-6 — 2) 
would require the disk to make up most of the rotation curve, which appears to be in conflict 
with the lack of faint stars or brown dwarfs in the solar neighbourhood (Gould, Bahcall, & 
Flynn 1996). As the effect of the disk is small, in this paper, we will simply adopt 

r disk ,- 6 « 0.35 - 1, (26) 

and completely ignore its z dependence. This range is indicated by two thick horizontal bars 
at the bottom left in Figures 3-6. 

3.3 Optical Depth of the COBE Bar: Red Clump Sources 

Having discussed the possible ranges of the bar parameters and the disk contribution, we 
can now investigate the optical depth of the COBE bar. For bright red clump giants in the 
bulge, we can apply eq. (|i~7l), which we rewrite as 
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(r (Qz))-21 M °- 5kpC ^MgMl^ [i | 1 (\ Z \Y]~^ (27) 
(T - 6(0 ' Z)) ~ 21 2xl O i°M ~ 0^64 0159 - [ + ^ U/ J ' 

where we have normalized the quantities £o(f3,n,k) and H(k,yo/xo,a) to the values for an 

oblate bulge, i.e., for yo/xo = 1 and n = k = 2. 

Table 1 gives the optical depths for four latitudes along the minor axis for the four models 

taken from Dwek et al. (1995). The field at (0°, —4°) is close to Baade window. The values 

are for M = 2.8 x 1O 1O M , at the upper end of eq. (|2~5"D, but this can be scaled linearly for 

other M. There are two things worth pointing out: first, at the same field the optical depths 

of different models differ by a factor 3-5; second, the optical depth drops with latitude for 

each model. From the Galactic center to Baade window, the decrease is about a factor of 

6 for the exponential diamond (ED) and the power-law ellipsoidal (PE) models, but only 

about a factor of 2-3 for other models. The fast decrease of the exponential diamond and 

power-law ellipsoidal models are obviously due to their steep density profiles (cf. Fig. lb and 

Table 1). 

These points are further illustrated in Fig. 3, where we have plotted the variation of the 
optical depth of the bar along the minor axis. The black dots with 2a errorbars are the 
inferred optical depths for the low latitude and high latitude clump giants (Alcock et al. 
1995b). Notice the data is suggestive of a decrease in the optical depth with latitude. 

Now let us examine Fig. 3 more closely. We notice that the errorbars are very large, which 
is hardly surprising because the two subsamples have only 10 and 3 events, respectively! 
Still, out of the four models, only the Gaussian ellipsoidal (GE) model can produce an 
optical depth as high as r_ 6 = 6. This is achieved only when the bar is very massive, with 
M = 2.8 x 1O 1O M and with a large disk contribution, 7di S k,-6 — 1- For the same condition, 
the oblate bulge and exponential diamond models are only at the 2a lower limit. The power- 
law ellipsoidal model produces an optical depth simply too small to be consistent with 
observation. This is mainly due to the large bar angle, a m 45°, rendering a very inefficient 
lensing configuration. 

3.4 Optical Depth of the COBE Bar: Main Sequence Sources 

We can carry out a similar analysis for the main sequence stars. The analytical formulae 
that apply to these stars are given in eqs. ( |G20| ) and ( |U21 ). But in Fig. 4, we only show 
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Figure 3. Optical depths on the z (minor) axis are shown for four COBE constrained bar models (columns 2 to 5 in Table 
1) for clump sources. The z coordinate and latitude are shown at the bottom and top axes respectively. The abbreviations 
for the models are given in Table 1. All the optical depths are for a very massive bar with M = 2.8 x 1O 1O M , but can be 
scaled linearly for other masses (the range is given in eq. ^fj). The approximate range of the optical depth contributed by the 
disk arc indicated with two short thick lines at the bottom left. The total optical can be calculated by first linearly scaling the 
plotted values for different M and then add the disk contribution. The two black dots are the optical depths observed for the 
low and high latitude clump giants (Alcock et al. 1995b). The 2tr errorbars are also plotted. The crosses indicate full numerical 
calculations for each model. 

the numerical results from eq. (0). Each shaded region corresponds to one model, with the 
upper and lower curves for (3 S — 1 and /3 S = 2 respectively. 

Quite strikingly the optical depth decreases with increasing (3 S) which may explain the 
smaller optical depth for main sequence stars f3 s > in the MACHO sample than for the 
clump giants j3 s = 0. For the Gaussian ellipsoidal model, the optical depth drops by a factor 
of 2 from f3 s = to (3 S = 2, while the power-law ellipsoidal model is almost unchanged. 
For the other two models, the decrease is about (20-50)%. Increasing (3 S effectively moves 
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Figure 4. Similar to Fig. 3, but for minor axis optical depth of four COBE models with source being bulge main sequence 
sources. Each shaded region corresponds to one model, with the top curve for j3 s = 1 and the lower curve for (3 S = 2. Again, 
a very massive bar with M = 2.8 X 1O 1O M is used. The optical depths for the full MACHO and OGLE samples are shown 
together with their 2a errorbars. 

the average source closer to us and reduces the optical depth by a factor proportional to 
f^-As; where X\ is the scale length in the line of sight (cf. eq.[|). The Gaussian ellipsoidal 
model has the largest major axis length and smallest angle a, hence the largest X\, the 
reduction for this model is therefore the strongest. The decrease of r for intrinsically faint 
sources can potentially provide interesting limits on the source luminosity function if the 
bar is sufficiently extended. 

In the same figure, the optical depths from 41 MACHO events and 9 OGLE events are 
shown with their 2a errorbars. With a bar mass of M = 2.8 x 1O 1O M and a disk contribution 
of 10 -6 , only the Gaussian ellipsoidal and exponential diamond models are within the 2a 
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errorbars of the MACHO sample. For the diamond exponential and power-law ellipsoidal 
models, the optical depths are too small even if we adopt a large bar mass and substantial 
disk contribution. 



3.5 Optimal Lensing Configuration of the Bar 

We have so far studied the COBE constrained bar models. Three out of the four models 
have difficulties in reproducing the observed high optimal depth. It is therefore interesting to 
examine whether these bar models - disregarding the COBE constraints- can produce the 
high optical depth at all. In other words, we want to know the maximum optical depth for 
each model by adopting the optical lensing parameters. For this exercise, we will concentrate 
on the clump giants as these are less prone to systematic effects such as blending and also 
because they appear to have larger optical depths. As the observed fields of MACHO and 
OGLE are centred near Baade window, we choose to maximize the optical depth at this 
window for each model. 

For each of the six sets of (/?, n, k) in Table 1, the optical depth at Baade window depends 
linearly on the total mass and on some combinations of the angle a, the axis ratio Uq/xq, and 
the vertical scale height zq. For the Baade window, the most favorable lensing configuration 
is when the vertical scale height zq = zbw ~ 590 pc (see the discussion at the end of §2.3.1), 
and the axis ratio Uo/xq = tana. 

In Fig. 5, we plot the optical depth for the above six bar models each at its most favorable 
configuration (zq = zbw an d Vo/ x o — tana) as a function of the bar angle a. Clearly for 
each model the optical depth drops quickly when a is increased from 0° to 20°, after which 
the curves become nearly flat. This is simply due to the factor 1/ sin(2a) in H(k, Uq/xq, a) 



(cf. eq. [18 



The six models in Table 1 have different radial profiles and boxyness due to their different 
(/3,n,k). Fig. 5 shows that the boxy bars with a flat density profile have an interesting 
signature of in optical depth. For any given bar angle the optical depth is the largest for 
the homogeneous rectangular bar. Intuitively this can be understood as follows. A large 
optical depth to the bulge requires putting as many lenses as possible at positions as distant 
as possible from the source. A model with a flat density profile in the line of sight places 
more lenses at the near side and more sources at the far side than a model with a steep 
density profile, hence a model with a flat profile is a more favorable lensing configuration. 
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10 20 30 40 

a(°) 



Figure 5. Optimal optical depths at Baade window are shown as a function of the bar angle a, for a very massive bar 
M = 2.8 X 1O 1O M0. To achieve the optimal optical depth, the vertical scale height is fixd to be zo = 590 pc, and the axis 
ratio yo/xQ = tana. The abbreviations for the models are listed in Table 1. The observed range of optical depth for the clump 
giants (Alcock et al. 1995b) is shown by two horizontal bars with arrows with the black dot indicating the central value. The 
approximate optical depth contributed by the disk lenses are plotted with two short thick lines at the bottom left. The total 
optical depth can be calculated by first scaling the plotted values linearly for different M and then adding the disk contribution. 

In addition, for a homogeneous bar the number of lenses does not decrease as one moves up 
on the minor axis before reaching the boundary. A rectangular bar furthermore stretches 
the distance between the lens and source even more as it has the longest contributing line 
of sight. A diamond shaped model, on the contrary, has the shortest contributing line of 
sight. These points become clear by comparing the length of segments of the dashed line 
intercepting the rectangular contour and the diamond contour in the bottom panel of Fig. 
la. Using the same argument, one can understand why the truncated boxy model (TB model 
in Fig. 5) also has a large optical depth. 
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This figure also points out the promise of the microlensing technique. If the present 
observed mean value of the optical depth is the true value, we see from Fig. 5 that all 
the bar models must have angle within 20° except the extreme homogeneous rectangular 
(HR) model. Unfortunately, if the lower 2a bound of the optical depth is true, and if one 
adopts a massive bar and a large disk contribution (r disk _ 6 = 1), then the current sample 
provides little constraint on the bar angle. The situation could be improved significantly 
with a factor of four increase in the sample size; the errorbar will be reduced by a factor of 
two. If the optical depth r_6 > 3.9 were confirmed, then the most plausible models would 
be a M > 2 x 1O 1O M truncated boxy (TB) bar with an angle a = 10° — 20°, and axis ratio 

< ^ ~ tan(a) plus a modest contribution of disk lenses. The optical depth r_6 ~ 12 
at the high 2a end of the MACHO clump giant sample is likely unphysical since it would 
argue for one of the following: (a) an extremely elongated bar with ^ = tan(a) « 0.05 and 
the line of sight almost coincides with the major axis of the bar: a = 2° — 5°; (b) an very 
elongated homogeneous rectangular bar with a ~ 10°; (c) a drastically different picture of 
the mass distribution in the Galaxy. The first two possibilities are probably ruled out by 
the COBE map, while the third one will have a host of difficulties with the gas and stellar 
kinematics. 

In Fig. 5, we have taken a as a complete free parameter, and the axis ratio is set such as 
to achieve the maximum optical depth. Now we want to examine possible constraint for the 
axis ratio. To do this, we take a = 10°, i.e., at the lower end of the inferred bar angle (cf. 
Table 1), keeping z = Zq, and study the variation of optical depth with the axis ratio. The 
results are plotted in Fig. 6. Big optical depth would imply that the axis ratio ^ is close to 
the optimal value tan(a). Unfortunately the current data gives only a rather weak limit at 
95% confidence level. However, if we take the mean value of the MACHO experiment and a 
disk contribution of 10~ 6 , we can place non-trivial constraints on the axis ratio. For example, 
for the GE model, we will limit the axis ratio to be 0.05 < yo/xo < 0.5; the observed value 
0.35 for GE model (cf. Table 1) is within this limit. More definite limits are only possible 
when the sample size is increased a few times and the errorbars are reduced. 

4 DISCUSSION 

We have evaluated the optical depths for a family of bar models. The analytical formulae 
developed are applicable to both the bulge clump giants and the main sequence stars. We 
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Figure 6. Optical depths at Baade window are shown as a function of the axis ratio ^ for a very massive bar, M = 

XQ 

2.8 X 1O 1O M , with a fixed angle a = 10° . The vertical scale height for each model is set to be zq = 590 pc. The abbreviations 
for the models are listed in Table 1. M = 2.8 X 10 10 Mq. The observed optical depth for all the clump giants by Alcock et al. 
(1995b) is shown by two horizontal bars with arrows with the central value indicated by a black dot. The approximate range 
of the optical depth contributed by the disk lenses are indicated by two short thick lines at the bottom left. The total optical 
can be calculated by first scaling the plotted optical depth linearly for different M and then adding the disk contribution. 

have studied the dependence of the optical depth on the bar parameters. With about 50 
microlensing events, we have demonstrated that microlensing already provides interesting 
constraints on the bar models. For example, the COBE constrained power-law ellipsoidal 
(PE) model produces an optical depth just too small to be consistent with the observed 
optical depths, primarily due to the large bar angle, a m 45°. The axisymmetric oblate bulge 
model is also only marginally compatible with the observation. The Gaussian ellipsoidal (GE) 
model clearly best reproduces the observation. It can produce r_6 = 6 if the Galactic bar 
is very massive (M = 2.8 x 1O 1O M ) and the disk contribution is Tdi S k,-6 = 1- We believe 
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that the truncated boxy (TB) model is also a viable model. It is therefore desirable to 
constrain the parameters for this kind of finite truncated bar using the COBE map. From 
our comparison of the optimal optical depths with the observations, we found the high 
optical depth already suggests a less than 20° angle between the long axis of the bar and our 
line of sight to the Galactic centre, and that the axis ratio of the bar is close to the optimal 
(cf. Fig. 5 and 6). 

In our comparison between the prediction and the observation, we have neglected a 
central nucleus and the dependence of r on the longitude. The longitude dependence in 
general reduces the optical depth. But as the mean longitude of the MACHO and OGLE 
fields are quite small, I < 2.5°, which corresponds to approximately 375 pc, the reduction is 
about 20% for the power-law ellipsoidal (PE) model and only about 10% for all the other 
models. In our Galactic bar models, the density profiles are flat near the Galactic centre, 
while a nucleus has long been observed in the infrared (Becklin & Neugebauer 1968). But 
since the MACHO and OGLE fields are all offset by more than 2° from the centre, a nucleus 
within 2° — 3° of the Galactic centre does not affect the microlensing prediction directly 
except that the nucleus contributes to the dynamical mass within 3 kpc. Including a nucleus 
would reduce the predicted optical depth outside the nucleus. But we estimate that the 
mass of the nucleus inside 2° is only about 5% of the total bar mass, much smaller than the 
uncertainty in the bar and disk mass. So the nucleus can be neglected for predicting optical 
depth at any field a few degrees 3°) away from the Galactic centre. 

In this study, we have used the MACHO and OGLE samples (Alcock et al. 1995b; Udalski 
et al. 1994). The number of analysed events (50) is still very limited, therefore the errorbars 
of the inferred optical depths are still quite large. The subsample of (13) clump giant events 
that we used has even larger errorbars. In addition, some of the long events in the MACHO 
sample are not well understood and these contribute about 1/3 of the optical depth (Han 
& Gould 1996). Therefore the current limits on the optical depth have to be interpreted 
with some caution. Fortunately, the number of events detected toward the Galactic bulge is 
increasing rapidly with time. In two years, the number of (analysed) microlensing events is 
likely to increase by four-fold. It will then become feasible to analyse the whole microlensing 
map of a large fraction of the bulge. If on the order of one hundred events are obtained at 
two or three low extinction fields close to the minor axis, e.g., the Sgr I field (1.4°, —2.6°) 
and Baade window (1°,— 3.9°), then we can measure the gradient for the optical depth, 
which we can use to distinguish bar models with different minor axis profile and boxyness 
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(cf. Fig. 3). If the optical depth remains high, then tighter limits on the bar parameters can 
be derived (cf. Figs. 5 and 6). By combining the microlensing map with the event duration 
distributions, it seems possible to disentangle the relative contributions of the disk and 
the bar. Furthermore, a comprehensive study combining the information in the gas, stellar 
kinematics, COBE map and microlensing surveys should provide precise values on all the 
bar parameters. We conclude that microlensing has become a promising and unique tool in 
probing the Galactic structures. 

This project is partly supported by the "Sonderforschungsbereich 375-95 fiir Astro- 
Teilchenphysik" der Deutschen Forschungsgemeinschaft. We thank Peter Schneider for com- 
ments on the paper. 



APPENDIX A: MORE DETAILS OF THE ANALYTICAL MODELS 



In this appendix, we give more technical details for the general family discussed in §2.3. 

For the density profile described by eqs. (|I5|) and ([LJ]), one can show that the central 
density is given by 
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where B is the complete Beta-function, and two approximations are given for K{k) and J(n) 
that are accurate within 3% for k > 1 and n > 1. 

The optical depth at any minor axis field depends on z as follows: 
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The optical depth toward the centre (r(0, 0)) and the optical depth for a source at the centre 
t(0, 0, 0) are given by 

(r(0, 0)) = r(0, 0, 0) 6(/3, n) = n)H(k, ^, a), (A6) 



c 2 2 
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where the quantities £q,£i,H are dimensionless functions of the bar angle, shape and radial 
profile. They are defined as 

£™dXf(X)f + °°dX 1 f(\X-X l \)X 1 
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where /(A) is defined in eq. ([16]), and the dependence on the bar angle and axis ratio has been 
collected in H(k, yo/xo, a). For most interesting combinations of (3 and n, £ x = (1 ± 0.1)\/2, 
and can be approximated within 5% by 
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£o also varies less than 30% from the value for a Gaussian elliposidal bar (£ = n 1//2 ) and 
can be approximated within 15% by 
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The Gaussian ellipsoidal model studied in §2.2 is a special case with (J3, n, k) = (oo, 2, 2). 
For this case, £j — > v2, £o = tt -1 ^ 2 - Another special case is ([3,n,k) = (0, oo,oo), which 
describes a rectangular homogeneous bar. The result for this model is given in Appendix B. 



APPENDIX B: RESULT FOR A HOMOGENEOUS RECTANGULAR BAR 

For a homogeneous rectangular bar, (/?, n, k) = (0, oo, oo), 

€i{P,n) = -, Co(/?,n) = -, H(k,—,a) = , (Bl) 

3 3 Xq 4 sin 2ct 



where w is defined in eq. (|A9|) . Towards the centre, 
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The optical depth at any field in the y — z plane 

(r(y, z)) = (r(0, 0)) , if \y\ < y ic and \z\ < z (B4) 

= (r(0,0))(^ yy^A ) if y ic < \y\ < y oc and \z\ < z 

\l-yic/yocJ 

= 0, otherwise, 

where 

yic = sin a — y cos a\, y oc = x sin a + y cos a, (B5) 

are the y coordinates of two of the four corners of the rectangular bar in the positive y side. 
Note (r(y, z)) is independent of y and z in a region near the centre. 



APPENDIX C: FIRST ORDER CORRECTION TO THE ANALYTICAL 
RESULTS 

Here we derive a first order correction of our analytical formula for the optical depth, taking 
into account of the finite distance to the bulge, and the source luminosity function. We show 
that the first-order correction is zero for the clump giants. 

The rigorous expression for the optical depth, as given by eqs. (1) and (2), can be 
rewritten as: 

(r} = ~ 2 J^n(D s )dD s ' dD. <x p{D a ) dD s , (CI) 

where the symbols are the same as in eqs. (1) and (2), and Ro is the distance to the Galactic 
centre. 

To the first order of (jjr^ > where a is a characteristic length scale of the bar, 
n / ,R + Xi (x s -Xi) 2 
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Changing the lens and source coordinates to the (x, y, z) frame, and assuming the density 
of source or lens is zero at the large radius, we have 

<r>»<n(»,«)> = *?£±£ (C4) 

where 

/+oo 
dx s p(x s ,y,z), (C5) 
-oo 
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h = 
and 

h = 2(1-/3.) 
= -2(1 -A) r 

J — < 

/ + GO >£ 
dx s — 1 [p(x s , y, z) - p(x s , -y, z) 
-OO lift 



+ 00 rj- 

dx s -^-p(x s ,y,z), 

-OO lhQ 

+ 00 <£ 

~b~p{ x s, -y,z), 



0, ify = 0orft 



dx. 



-OO 
OO 



dx,, 



[I- (3 S 



(3 S 



-OO 
X s 

-OO 
OO 

-OO 
OO 



dXl P+y^+y 
dX\ P—y^—y 



X s Xl)_ + 2 ( 1 _ /3s ^s{X s Xl 



Rq 

(x s - xi) 2 



+ 2(1 -A) 



-Ro 

Xl(Xj - Xg) 
Rq 



dx s 
dx 



dx\ p dd 

s 



/y»2 /y»2 



2i? 

(x s - Xi) 2 



(C6) 

(C7) 

(C8) 

(C9) 
(CIO) 

(Gil) 

(C12) 

(C13) 



-A 



(X s ~ Xi) 
Rn 



if 2/ = 0, 



2i? 

^ p(x s ,0,2;) p(z z ,0, z) 

J-oo J —OO 

where 

P+2/,+2/ = P( x ^ y, z) p(^, y, z), P-y-y = p(x s , -y, z) p(x h -y, z), 

Peven = P+y,+y + P-y-y, Podd = P+y,+y ~ P-y-y 

In the above we have used (ri(y, z)) to denote the first order corrected optical depth. If 
we use {rn(y, z)) to denote the optical depth given by eq. (f|) in the main text, then (ri(y, z)) 
reduces to (ro(y, z)) if the first order correction terms I± and I3 are set to zero. 

Eqs. ( |C8| ) and ( |U 1 2| ) are derived by applying a change of the dummy variables 



(C14) 

(C15) 
(C16) 



x s -> -xi, xi -> -x s , (C17) 
and the central symmetry of bars 

p(-x,y,z) = p(x, -y,z), (C18) 

to eqs. (]C7| ) and ( |(J11| ). We then obtain eqs. ( |C5| ) and ( |C13| ) after adding eqs. ( |CTD and 
(p8D, and eqs. ( |Cil| ) and ( |<J12| ), and then dividing both sides by two. Eq. QU14D is because 
Podd = on the minor axis. 

A nice result from this is that for red clump stars with presumably j3 s = 0, I± = I3 = 
for any line-of-sight along the minor axis y = 0. Hence (ti(0,z)) = (to(0,z)), namely, the 
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optical depth given in the main text is correct to the first order. This conclusion is true for 
any bar density profile, shape and orientation. 

For main sequence sources with f3 s > 0, our formula over-estimates the optical depth 
and has to be corrected down. For the density profiles described by eqs. (|T5| ) and (|T6|), after 
some tedious algebra, the optical depth, corrected to first order, is given by 

fp{0,0,z)Y 



(n(o,z)} = (t (o,*))(i + ^) = <toM) 



1 - UPs 



Ro Vp(o,o,o). 



(7-0(0,*)) 



Ra 



1 + 



/3-3 



where 



cos fc a 



k 

;r 



+ 



sm k a 



(C19) 
(C20) 

(C21) 



The expression for ^(f3,n) ~ 3« is accurate within 15% for most bar models, and is exact 
for homogeneous {(3 = 0) bar. For Gaussian ellipsoidal bars, £4 = \fTi. Other exact values for 
£4 for a few models are given in Table 1. From eq. (|U20p, the first-order correction is bigger 



away from the plane for power-law models {(3 > 3), smaller for truncated bars {f3 < 0) and 
independent of z for (3 — > 00 (e.g., the Gaussian and exponential models). 

We have also checked these formulae with fully numerical calculations. We find that the 
above first order corrected results hold within 5%-15% for the full range of models discussed 
in the main text, which have a/Ro < |. 

For readers' convenience, we summarize in the following how to use our analytical for- 
mulae: 

(i) Find the value of t;o(j3,n). £0 for six comibinations of (/3, n) are tabulated in Table 1. 
For combinations not found in the table, the approximation given by eq. (|A11|) should be 
sufficient for most cases. When more accuracy is desired, use eqs. (|A7|) and (|A8|) . 

(ii) Use eq. ([A|]) to calculate the function H(k,yo/x ,a). 

(iii) Now use eq. (|27|) to calculate the zeroth-order optical depth. 

(iv) The first order corrected optical depth can be obtained using eqs. (|C21|) and (|C20|). 
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